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Recent numerical and analytical work has shown that for the square-lattice Heisenberg model 
the boundary can induce Dimer correlations near the edge which are absent in spin-wave theories 
and non-linear sigma model approaches. Here, we calculate the nearest-neighbor spin correlations 
parallel and perpendicular to the boundary in a semi-infinite system for two different square-lattice 
Heisenberg models: (i) A frustrated Ji — J2 model with nearest and second neighbor couplings and 
(ii) a spatially anisotropic Heisenberg model, with nearest-neighbor couplings J perpendicular to 
the boundary and J' parallel to the boundary. We find that in the latter model, as J' / J is reduced 
from unity the Dimer correlations near the edge become longer ranged. In contrast, in the frustrated 
model, with increasing J2 , dimer correlations are strengthened near the boundary but they decrease 
rapidly with distance. These results imply that deep inside the Neel phase of the Ji — J2 Heisenberg 
model, dimer correlations remain short-ranged. Hence, if there is a direct transition between the 
two it is either first order or there is a very narrow critical region. 

PACS numbers: 



I. INTRODUCTION 

Square-lattice antiferromagnets have been studied ex- 
tensively in recent yearsi^^&L&ii. Yet, new surprises 
continue to arise. In particular, recent Quantum Monte 
Carlo studies^ by Hoglund and Sandvik have shown that 
the existence of a free edge induces pronounced Dimer- 
ized correlations in the system. In a follow up workii it 
was shown by Metlitski and Sachdev that the presence 
of a boundary induces dimer correlations perpendicular 
to the boundary. And since the correlations decay with 
distance from the boundary, their gradient induces alter- 
nation in the spin-correlations parallel to the boundary, 
leading to specific pattern of nearest-neighbor spin cor- 
relations observed by Hoglund and Sandvik in their sim- 
ulations. These effects are absent in spin-wave theories 
and non- linear sigma model approaches. 

Over the past few years there has been considerable in- 
terest in the possibility of direct continuous phase tran- 
sitions between Neel and Valence Bond Crystal (VBC) 
phasesi^^^^iii^^^i^i. Such phase transitions 
have been called deconfmed quantum criticality and are 
marked by the liberation of spin-half degrees of freedom 
as well as the existence of massless spin-singlet photon 
field. Strong numerical evidence for such a scenario has 
been provided in Sandvik's J-Q mode l 22 ' 23 ' 24 , where the 
Heisenberg model is supplemented by a 4-spin interac- 
tion around a plaquettc. An alternative possibility of a 
weakly first order transition has also been raise d 25 ' 26 . 

A more realistic model of two-dimensional square- 
lattice quantum antiferromagnets is the spin-half J± — J2 
model, where there is nearest-neighbor interaction J\ and 
second neighbor interaction J2. Increasing J2 increases 
spin frustration and is known to lead to a magnetically 
disordered state at intermediate J2/J1 values 5 -^. There 
is substantial and growing body of numerical evidence 
that the magnetically disordered phase has Valence Bond 



Crystal (VBC) order^i^. The question of whether the 
transition between the Neel and Dimer orders is contin- 
uous or first order remains a subject of debat o 27 ' 28 . 

Here, we would like to use the edge induced dimer 
correlations as a probe of growth of dimer correlations 
inside the Neel phase and thus address the possibil- 
ity of a diverging dimer correlation length in the Neel 
phase. We study two models. A spatially anisotropic 
model with interactions J and J' along the two axes. 
We choose the boundary to be parallel to the direction 
of the weaker coupling J'. It is well known that one- 
dimensional Heisenberg model has power-law decaying 
Valence Bond correlations. Thus as one approaches the 
limit of small J' one expects to see the edge induced cor- 
relations to have a long length scale. This model acts as 
a test case for our method. We also study the J\ — J2 
Heisenberg model. It is for this model that one would 
like to see how the range of dimer correlations grows near 
the boundary as spin-frustration given by the parameter 
J2/J1 increases and one approaches the phase transition, 
where Neel order is lost. 



II. SERIES EXPANSION 

The antiferromagnetic Heisenberg models defined by 
two coupling constants Ji and Ji (or by J and J') are 
shown in Fig. [T] We consider a semi- infinite system, with 
a boundary parallel to the V-axis also shown in Fig. 1. 
Since we are considering a system inside a colinear Neel 
ordered phase, we develop an Ising series expansion 2 ^, 
where all Heisenberg couplings are written as 



QZ QZ , \ ( qx qx i qV QV \ 

Oi Oj ~r Ayo^ Oj ~f Oj j. 



(1) 



The parameter A acts as an expansion parameter. We 
develop series expansions for on-site local magnetiza- 
tion (Sf) as well as for nearest neighbor spin-correlations 
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FIG. 1: (left) J\ — J 1 model on the semi-infinite square lat- 
tice. The interaction Ji = 1 is perpendicular to the edge, 
the interaction J' is parallel to the edge and it is allowed to 
vary in the range < J' < 1. (Right) Ji — J2 model on the 
semi-infinite square lattice. J\ = 1 and < J2 < 0.3. The 
parameter R defines the distance of a particular lattice plane 
from the edge. 



(S; • Sj), parallel and perpendicular to the boundary. In 
the semi-infinite system, these quantities depend on the 
distance R from the boundary. In the series expansion 
method, the boundary can be accomodated by account- 
ing for the graphs that terminate at the boundary. Apart 
from this, the formalism of linked cluster expansions re- 
mains unchanged. 



bulk value. In the non-linear a model and spin wave the- 
ory the change in magnetization follows a l/R behavior. 
On general grounds one expects the non-linear a model 
results, when expressed in terms of renormalized param- 
eters, to be exactii for large-R. The reduction is less in 
our calculation upto the largest distance studied, that is, 
R = 5. Part of the reason maybe that the asymptotic 
behavior may set in at significantly large-R due to the 
dimer-correlations at the boundary. However, it is also 
likely that the uncertainty in our calculations are much 
larger than shown. Our estimate of the bulk magneti- 
zation is 0.302. If we replace it by the more accurate 
results from higher order series expansions^ or quantum 
Monte Carlo simulations 22 , which is 0.307 it would shift 
our calculated curves up by 0.005 and bring them closer 
to the spin-wave results. This discrepency in the bulk 
estimates implies that the uncertainties are much bigger 
than estimated by the fits and they are particularly mag- 
nified at larger R because we are taking the difference of 
two quantities which are close in magnitude. 

In Fig. [51 we show the nearest-neighbor spin- 
correlations perpendicular to the boundary for the J — J' 
model. This is the direction of the stronger coupling. In 
the one-d limit, one expects the free end to induce dimer 
correlations in the system that decay as a power-law away 
from the boundary. Indeed, we find that as the system 
becomes more and more anisotropic, the dimer correla- 
tions become more and more long ranged. 



III. RESULTS: CORRELATIONS AND 
EXCITATIONS NEAR THE EDGE 

First we present the results for the square-lattice 
Heisenbcrg model. The nearest neighbor correlations 
parallel and perpendicular to the boundary are shown 
in Fig. [21 These are obtained by d-log Pade approximant 
analysis of the series. They agree well with the results 
of Sandvik and HoglundiS. The important thing to note 
is that they both decrease rapidly with distance and by 
R = 5 they differ from the bulk value by less than 0.1%. 

The on-site magnetization also changes near the 
boundary. The results for magnetization are more sen- 
sitive to extrapolation methods than spin-spin correla- 
tions because one expects a square root singularity for 
this quantity. This means that contributions of higher 
order terms only decay as 1/y/N. We have followed the 
method used in Ref. [3fJ for the bulk system. We obtain 
partial sums Sn of series coefficients and then fit them vs 



in Fig. [3] for values of R > 2. We deduce the uncertainty 
in the magnetization by the uncertainty in the linear fits. 

Results obtained this way are plotted in Fig. [4] where 
they are compared to the non-linear a model and spin 
wave results^. The on-site sublattice magnetization is 
diminished at the edge and its reduction is comparable 
to what is obtained in spin-wave theory. Away from the 
edge the sublattice magnetization should approach its 
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FIG. 2: (Color online) Correlation function of the spin- 7; 
Heisenberg model on the semi-infinite (SI) square lattice for 
bonds parallel and perpendicular to the edge as a function 
of distance R. Notice that bonds perpendicular to the edge 
have their centers at half-integer values of R. The y-axis is 
— Sl ' >SI — s j)°° where the 00 index refers 



to the bulk value. 



In Fig. [6j we show the nearest-neighbor spin- 
correlations for the Ji — J2 model. Correlations both 
parallel and perpendicular to the boundary are shown. In 
this case, we find that while a frustrating second neighbor 
interaction enhances the dimerization near the boundary, 
it does not appear to increase the range over which dimer 
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FIG. 3: (Color online). Partial sums of series expansion co- 
efficients for the on-site sublattice magnetization of the semi- 
infinte square lattice model. The fit for different values of the 
parameter R are shown. See text for details 
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FIG. 4: On-site sublattice magnetization for the spin-^ 
Heisenberg model on the semi-infinite square lattice. The 
non-linear sigma model and spin wave results from Ref. [ll| 
are also shown. The i/-axis is AM = (Msi — Moo). 



correlations extend. The convergence of our analysis be- 
comes poor as we get close to the bulk transition away 
from Neel order, which has been estimated to be in the 
range J 2 /Ji « 0.35 - OA^l 

These results show that in the J\ — J2 square-lattice 
Heisenberg model, one does not have appreciable range 
Valence Bond Correlation in the bulk even with signif- 
icant frustration. They suggest that a direct transition 
between Neel and Dimer phases is likely first order. Our 
study can not rule out the possibility that the dimer cor- 
relations build up very quickly close to the transition. 
This would imply a very narrow critical region in this 
model. 

We have also calculated the spin-wave spectrum for 
the magnon states that are bound to the surface for the 
nearest-neighbor square-lattice Heisenberg model. The 
momentum parallel to the surface is a good quantum 
number. In the series expansion calculations, the spin- 
flip states right at the boundary have a different excita- 
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FIG. 5: (Color online) Correlation function of the spin-i 
Ji — J' model on the semi-infinite square lattice for bonds 
perpendicular to the edge as a function of distance R for se- 
lected values of J'. The quantity ACVj/Cfj shown on the 
y-axis is defined in the caption of Fig [2] 
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FIG. 6: (Color online) Correlation function of the spin-| Ji — 
J2 model on the semi-infinite square lattice for bonds parallel 
(top panel) and perpendicular (bottom panel) to the edge as a 
function of distance R for selected values of J2. The quantity 
ACij/Cfj shown on the y-axis is defined in the caption of 
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FIG. 7: Excitation spectrum of the Heisenberg model on the 
semi-infinite square lattice for < k x < rr, where k x is the 
momentum parallel to the edge. The bound states found by 
spin- wave theory and series expansion calculations are shown. 
The thin solid line represents the upper limit of the continuum 
of states. 

tion energy from those which are away from the bound- 
ary. Thus these states get separated from the bulk states 
starting in zeroth order. Upon extrapolation to the 
Heisenberg model, we find the dispersion of these sur- 
face magnons as shown in the Figure [7J Also, shown are 
results from the spin-wave calculations of Metlitski and 
Sachde\*ii. The latter has been renormalized to have the 
same spin-wave velocity as the bulk. Our results are in 



agreement with the latter that for a large part of the 
Brillouin zone, the surface states are hugging the contin- 
uum. Only very near k — 7r/2, they clearly separate from 
the continuum. In this region, the binding energy in our 
calculation is smaller than in spin-wave theory. 



IV. CONCLUSIONS 

In this paper we have studied the spin-correlations and 
excitations near the boundary of two dimensional Heisen- 
berg antifcrromagnets. Two different square-lattice mod- 
els are considered. One where the exchange coupling par- 
allel to the boundary is smaller than those perpendicular 
to the boundary. In this model, we find that the bound- 
ary induced dimerization becomes more and more long- 
ranged as the anisotropy is increased. The second model 
is the Ji — J2 Heisenberg model, with nearest and sec- 
ond neighbor exchange interactions. In this case, we find 
that as frustration is increased in the model, the bound- 
ary induced dimerization increases close to the boundary 
but its range does not change significantly. This suggests 
that in the J\ — J2 model, the Neel phase does not de- 
velop long-range dimer correlations. Hence, either the 
transition from Neel to dimer order is first order in this 
model, or there is a very narrow critical region. 
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